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TECHNICAL NOTE 2868 


REFLECTION OF A WEAK SHOCK WAVE FROM A BOUNDARY 


LAYER ALONG A FLAT PLATE 


I - INTERACTION OF WEAK SHOCK WAVES 
WITH LAMINAR AND TURBULENT BOUNDARY LAYERS 
ANALYZED BY MOMENTUM-INTEGRAL METHOD 


By Alfred Ritter and Yung-Huai Kuo 
SUMMARY 


The present paper is concerned with the phenomena encountered when 
a plane oblique shock wave is incident upon the boundary layer of a flat 
plate. In an effort to simplify the problem, the flow field was divided 
into a viscous layer near the wall and a supersonic potential outer flow. 
The pressure disturbances due to the incident wave would be propagated 
upstream and downstream in the subsonic portion of the boundary layer, 
thus giving rise to perturbations of the boundary layer. By restricting 
the study to infinitesimal incident compression waves, only small per- 
turbations were encountered and hence the ordinary linearized theory 
could be applied to the outer flow. In the laminar case, the boundary- 
layer treatment was based upon a momentum-integral equation previously 
derived by Howarth. The two flows must be compatible; hence, the deflec- 
tion of the streamlines near the boundary layer was expressed in terms 
of the vertical velocity component along the edge of the boundary layer 
and this relation was used as a boundary condition for the outer flow. 
The boundary condition determined the form of solution upstream and down- 
stream of the point of incidence. Determination of the constants of 
integration was accomplished by a consideration of conditions at infinity 
and a matching of the two flows at the point of incidence. With the 
outer flow thus determined, boundary-layer growth and pressure distribu- 
tion were computed and results for the laminar case were obtained as 
follows: 


i (a) The pressure disturbance along the wall decreased exponentially 
from a definite value at the point of incidence to zero far upstream of 
the point of incidence. Downstream of the point of incidence, the pres- 
sure rose to a maximum value and then dropped off to the value corre- 
sponding to regular reflection. 
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(b) The disturbances produced by the interaction decayed exponen- 
tially upstream; for a free-stream Mach mmber of approximately 2 and a 
Reynolds mmber of approximately 1500 in the undisturbed boundary-layer 
displacement thickness the upstream influence was of the order of 30 
boundary-layer displacement thicknesses. ; 


(c) The "self-induced" pressure gradient along the wall was such 
that the boundary layer might separate ahead of the point of incidence. 
If separation occurred, the separation point moved upstream as the shock 
strength was increased. With increasing Reynolds number, the separation 
point also moved upstream, whereas for increasing Mach number, the sepa- 
ration point moved downstream. 


In the turbulent case the upstream influence was quite small and the 
incident wave must be reflected as a shock wave. 


INTRODUCTION 


It has been found that if the free-stream subsonic Mach number 
becomes high enough so that local supersonic zones are formed on an air- 
foil, sharp changes in the airfoil characteristics occur which cannot be 
explained by classical aerodynamics. The resulting loss of lift is 
accompanied by a large increase of drag in consequence of the appearance 
of shock waves on the surface of the airfoil. However, close study 
reveals that the drag increase is too large to be accounted for by the 
shock loss and change in skin friction. It must be caused by the sudden 
change of the flow pattern. This seems to indicate that the shock wave, 
when formed over the airfoil surface, modifies the character of the 
boundary layer in such a way as to create a wider wake. 


At Guidonia, Ferri, by examining the measured pressure distributions 
over airfoils at supersonic speeds, found that for the forward portions 
of the airfoil the experimental pressure distribution agrees quite well 
with that calculated from potential-flow theory (reference 1). As the 
trailing edge is approached, however, the experimental pressures along 
the upper surface, for positive angles of attack, become considerably 
higher than the calculated ones. (At negative angles of attack, this 
behavior occurs along the Mower surface.) This discrepancy is under- 
standable in the light of the fact that the boundary layer on the rear 
portion of the airfoil is extremely sensitive to pressure disturbances. 
Since the flow far away from the airfoil has to return to its original 
pressure and direction, a shock wave must emanate at or near the trailing 
edge. There is a sharp pressure rise across the shock and pressure dis- 
turbances are transmitted upstream in the subsonic portion of the boundary 
layer causing a thickening of the layer near the trailing edge. This 
thickening, in turn, generates compression waves which travel downstream 
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and interact with the shock wave. The net result is that there is, 
starting from some point forward of the trailing edge, a gradual compres- 
sion to the main-stream pressure. The pressure distribution is hence 
altered in such a manner as to give values of lift and pressure drag 
smaller than those calculated from theory. The measured moment coeffi- 
cient will also differ from the theoretical values. 


The preceding discussion shows that, in the transonic and supersonic 
regimes, the Mach number alone is insufficient to determine the flow 
Characteristics; the Reynolds number can also be important. This was 
first demonstrated by Ackeret, Feldmann, and Rott in,experiments which 
established the close relationship between the shock-wave pattern and 
the Reynolds number (reference 2). These prove beyond doubt that the 
flow far away from the wall depends intimately upon the character of the 
boundary layer, that is, the Reynolds number. Thus, under such circum- 
Stances, the concept of the boundary layer requires modification. 


As a first step, a simplified problem of a plane shock wave incident 
upon the laminar boundary layer over a flat surface will be dealt with. 
For a number of years, Dr. H. W. Liepmann and his associates at the 
California Institute of Technology have been conducting experimental 
studies of shock-wave boundary-layer interaction. Recently they investi- 
gated the problem of the reflection of shock waves from boundary layers 
(reference 3). Based on the gualitative experimental data, some impor- 
tant conclusions were drawn regarding the reflection of &n incident shock 
wave from & boundary layer: 


(&) The type of boundary layer, whether laminar or turbulent, 
markedly &ffects the interaction. With a turbulent boundary layer, the 
reflection is practically the same as that for a regular reflection in 
nonviscous flow theory. For a laminar boundary layer, however, there is 
a large interaction zone near the point of incidence that is guite dif- 
ferent from the regular inviscid reflection. The incident wave is appar- 
ently reflected as from a constant-pressure surface in the form of a 
Prandtl-Meyer fan. In returning to its final direction parallel to the 
wall, the deflected flow, behind the expansion, recompresses to the pres- 
sure appropriate to that behind a regular reflection. 


(b) For M = 1.4 and Rex = = = 0.9 x 106, the influence of the 


incident wave extende upstream in a laminar layer a distance of the order 
of 50 boundary-layer displacement thicknesses. On the other hand, the 
upstream influence is practically negligible for the turbulent boundary 
layer. 


(c) It appears that, except for very weak incident waves, the laminar 
layer almost always separates somewhere in front of the point of inci- 
dence, whereas no separation was found for the turbulent layer. 
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Quite recently, Barry, Shapiro, and Neumann (reference 4) made a 
similar study and their results are in substantial agreement with the 
GAICIT results. An advantage of their experiments, however, is that, 
although the tests were carried out at constant Mach number, the Reynolds 
number and shock strength were varied so that more qualitative as well as 
quantitative data are available to serve as & guide to future theoretical 
work. 


Theoretical solutions of this problem have been almost as meager as 
the experimental results. A first attempt was made by Howarth (refer- 
ence 5)5 who considered the case of & wave incident upon the interface 
bounding two semi-infinite uniform streams, one supersonic and one sub- 
sonic. Viscosity and heat conduction were neglected and the equations 
were linearized. He then demonstrated the upstream propagation of dis- 
turbances in the subsonic portion of the flow field and showed that an 
incident compression wave is reflected as compression upstream of the 
point of incidence while being reflected as expansion downstream of the 
point of incidence. However, these results can only be regarded as quali- 
tative because of the obvious shortcomings of the model. In an effort 
to make the Howarth model more realistic, Tsien and Finston simulated 
the boundary layer with & uniform subsonic stream bounded by & wall on 
one Side and a semi-infinite, uniform, supersonic stream on the other 
(reference 6). As before, viscosity and heat conduction were neglected 
and only small disturbances were considered. For the case of a compres- 
sion wave incident upon the interface separating the two streams, it was 
shown that except for the local interaction, the incident compression 
wave is regularly reflected. Locally, however, the pressure along the 
interface exhibits compression ahead of the point of incidence and expan- 
sion immediately behind it. This local condition is qualitatively the 
same as that observed in experiments with shock reflection from laminar 
boundary layers; consequently there arose the speculation that perhaps 
the effects of viscosity and heat conduction were actually not too impor- 
tant in comparison with the effect of coexistence of supersonic and sub- 
sonic streams. If this were the case, it should follow that an improve- 
ment in the Tsien-Finston model could give rise to more quantitative 
results. Now an obvious discrepancy between the Tsien-Finston model and 
the physical situation is that the uniform subsonic stream, in attempting 
to simulate the boundary layer, is incapable of satisfying the "no-slip" 
condition at the wall. Thus, the next logical step in the development 
of the theory would appear to be that of considering a main-stream velo- 
city that varies from zero at the wall to a given uniform supersonic 
velocity a short distance from the wall, and, in fact, this has been con- 
sidered by Robinson (reference 7) and Lighthill (reference 8). Robinson 
assumed that the main-stream velocity varies continuously from Zero at 
the wall to supersonic velocity some distance from the wall, whereas 
Lighthill assumed that the Mach number varies from zero at the wall to 
some supersonic value a short distance from the wall. A weak wave is 
incident upon the boundary layer and the reflected waves and upstream 
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influence were evaluated. Now with the improved model, however, both 
investigators found, contrary to experimental evidence, that the upstream 
influence is negligible. Although Robinson was chiefly concerned with 
the determination of the upstream influence, Lighthill considered, in 
some detail, the local reflection of the shock wave. He concluded that 

a shock is reflected locally as a "pressure ridge," that is, a rapid 
compression followed immediately by a rapid expansion. Now a note was 
added in the proof of reference 8 to the effect that experimental results 
of Mair and Bardsley (reference 9), obtained at the Fluid Motion Labora- 
tory, Manchester, establish that the conclusions of negligible upstream 
influence and the shock being locally reflected as a pressure ridge are 
correct for the reflection of a weak shock from a turbulent boundary 
layer. (These results, however, had been previously observed by Ackeret 
and Liepmann.) As Lighthill pointed out, the conjecture that the theory 
is correct for velocity profiles typical of turbulent boundary layers 
but incorrect for profiles typical of laminar boundary layers may be 
valid. But until such time as the theory is modified to account for vis- 
cosity, as suggested by Lighthill, the present theory is incapable of 
predicting the effects of the interaction between a shock wave and a 
laminar boundary layer. 


The preceding theories have neglected viscosity and heat conduction 
while considering infinitesimal waves and small disturbances. It would 
appear that the assumption of infinitesimal waves and small disturbances 
is valid since the linear theory, to be a good approximation, requires 
that the slopes of the streamlines be small. For weak incident shock 
waves, the slope of the displacement thickness would probably be small, 
except perhaps in the immediate vicinity of the wave. The neglection of 
viscosity, on the other hand, seems to be more serious. If the effects 
of viscosity were not too important, one would expect that, for fixed 
Mach number and shock strength, the effect of Reynolds number should be 
rather small. Actually, this is not the case since the results of Barry, 
Neumann, and Shapiro clearly show that changes in Reynolds number have a 
marked effect upon the shock-wave - laminar-boundary-layer interaction 
(reference 4). Hence it appears that & theory capable of predicting the 
effects of this complicated phenomenon must include the effects of 
viscosity. 


At present, as far as the authors are aware, the only solution that 
includes viscous effects was given by Lees (reference 10). The Lees 
theory is based on an approximate relationship between the slope of the 
displacement thickness and the external pressure gradient derived from 
Prandtl-Meyer flow. By combining this relation with the momentum- 
integral equation and neglecting terms of the order of Rerl Lees arrived 
at a third-order, linear, differential equation for the pressure. On the 
assumption that this equation is valid all along the boundary layer, he 
was able to determine the boundary-layer growth and pressure distribution. 
He found the upstream influence comparable with that which is observed 
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experimentally and showed, also, that the boundary layer always sepa- 
rates except in those cases where the incident wave is quite weak. How- 
ever, his theoretical pressure distributions fail to exhibit the charac- 
teristic downstream behavior. This point will be discussed later. 


Following the approach of Oswatitch and Wieghardt in their study of 
the growth of disturbances in a supersonic stream outside laminar or 
turbulent boundary layers (reference 11), the problem is treated as an 
outer flow, with a shock, in equilibrium with a boundary-layer flow. 
Instead of expressing the deflection of the streamline in terms of pres- 
sure rise, as in Lees' problem, a procedure after Oswatitch and Wieghardt 
is taken by connecting the vertical velocity of the outer flow with the 
streamline deflection, as given by the momentum integral of the boundary- 
layer flow. This reduces the problem to an inviscid one which can be 
solved systematically to any order of approximation. 


This study was conducted at the Graduate School of Aeronautical 
Engineering of Cornell University under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 
a speed of sound 
A,B,C,D,E constants 
b = 51*/5,* 
Cp) Cy specific heats at constant pressure and constant volume, 
respectively 
f.g functions of E - mo] and E + mn, respectively 
l 


F = 27 - 2 +" 


81» 80.8398) defined by equations (8) 
I 
G = zw - DE 


h function of E - mam 


H =6*/0 
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coefficient of thermal conductivity 


Karman momentum equation for compressible viscous fluids 


Mach number (U/a) 


pressure 
perturbed pressure 

Prandtl number (ucp/ x) 
universal gas constant 


Reynolds number (U50*/v0) 


temperature 


velocity components parallel and normal to flow direction, 
respectively 


velocity perturbations nondimensionalized with U 
free-stream velocity 


coordinates parallel and normal to flow direction, 
respectively 
values of x and y in transformed: plane [o = Xj 
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constant 


2 
a — 
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y ratio of specific heats (cp/cv) 

5 boundary-layer thickness 

O*. boundary-layer displacement thickness 

Og* undisturbed boundary-layer displacement thickness 

O' boundary-layer thickness in transformed plane 

§'* boundary-layer displacement thickness in transformed plane 
A disturbed boundary-layer displacement thickness nondimen- 


sionalized with Sof 
€ flow-deflection angle 


n = y/5' or y/50' 


9 boundary-layer momentum thickness | H i (2 ~ = pu x) 
| O Ue} Pele 

g' boundary-layer momentum thickness in transformed plane 

A1 A23 X3 roots of equation (16) 

Às SL iy kedā. ne?) e 

Yo 2 dx 

" coefficient of viscosity 

V kinematic coefficient of viscosity 

E = x/o,* 

p density 

sal LE (19232 42) 

T Shear stress 

0 perturbation velocity potential, nondimensionali zed 


with Udo” 


2 
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o total potential 

X function of x' and y 

y stream function 

Subscripts 

d evaluated at distance of upstream influence 

e evaluated along edge of boundary layer 

1 l2 2, OE 3 

L laminar case 

O standard condition, such as stagnation point or undis- 
turbed state 

8 Separation point 

t turbulent case 

W evaluated at wall 

136,3 evaluateā in regions 1, 2, or 3, respectively (see fig. 1) 

co evaluated in undisturbed free stream 


OUTLINE OF PRESENT INVESTIGATION 


In the present study the effect of an infinitesimal compression 
wave incident upon the laminar boundary layer along a flat plate in 
supersonic flow will be considered. The pressure gradient in the flow 
direction is determined to the order of approximation of the ordinary 
boundary-layer theory solely by the shape of the boundary. The poten- 
tial flow thus completely determines the boundary-layer growth. In the 
problem of interaction between & shock wave and & boundary layer, how- 
ever, the phenomena are quite different. In that case the pressure dis- 
turbances that are propagated upstream in the subsonic portion of the 
boundary layer will affect the boundary-layer growth. It becomes clear 
that the outer flow influences the boundary-l&yer growth and that the 
boundary-layer growth, in turn, influences the outer flow, so that a 
solution that is simultaneously compatible with the two flows must be 
sought. E 
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Another significant difference between the interaction and noninter- 
action problems occurs in the separation phenomenon. For the usual case 
of boundary-layer flow against a pressure gradient, it is known that the 


flow separates from the wall at the point where 3 = O. Downstream 
W 
of the separation point the original flow is deflected away from the wall 
by the backflow, a vortex layer separating the two flows. The interac- 
tion of a shock wave with a laminar boundary layer, however, produces a 
different form of separation. Experimental evidence indicates that along 
the wall there exists a short thin region, extending slightly upstream 
from the foot of the shock wave, in which the flow has practically zero 
velocity; the boundary between this separated flow and the main boundary- 
layer flow is a vortex sheet. A plausible explanation of this phenomenon 
is proposed by Lighthill, who suggests that the pressure discontinuity 
occurring when a shock interacts with a laminar boundary layer causes the 
flow to separate into some form of "bubble" at the base of the shock wave. 
Once such'a bubble is produced, the boundary layer upstream of the sepa- 
rated region will be deflected so as to increase the external pressure 
gradient. This results in further separation of the boundary layer and 
hence further upstream deflection. This process of repeated separation 
causes the edge of the bubble to move upstream until such time as the 
induced pressure gradient is no longer able to cause separation of the 
boundary layer. The net result is the one indicated in experiments. 
Immediately behind the point of interaction, the deflected flow, in 
returning to its original direction parallel to the wall, undergoes a 
very rapid compression, thereby causing an extremely large pressure gra- 
dient. It is known from experimental results that the boundary layer 
downstream is turbulent, but whether the extreme pressure gradient causes 
separation and hence causes the boundary layer to become turbulent or 
whether the flow at that point is already turbulent and thus can sustain 
large pressure gradients must remain at this time a matter of conjecture. 
Whatever the case may be, in the problem for weak shock the flow will be 
assumed to remain laminar downstream of the point of incidence (cf. ref- 
erence 12). Positive infinity in relation to the problem will be taken 
to mean a distance downstream of the order of 100 boundary-layer thick- 
nesses. In addition, any mention of boundary-layer separation will refer 
to the "dead air bubble" phenomenon. 


In the present theory it is assumed that the boundary-layer equa- 
tions are applicable over the whole plane. Strictly speaking, the valid- 
ity of the usual boundary-layer assumptions near the base of the shock 
wave is doubted, since separation occurs. However, it is believed that 
for weak incident waves, a theory based on the above assumptions will 
still enable one to determine some of the important characteristics of 
the flow. In particular it is desired to determine the extent of the 
separated region, since this investigation is merely the first step of a 
proposed theory which will endeavor to take separation into account. 
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The more refined theory is then expected to furnish, in the separation 
zone, slight corrections to the simplified theory. 


The flow field is assumed to be divided into two distinct regions 
as follows: 


(a) A thin layer adjacent to the plate wherein there is laminar- 
type boundary-layer flow 


(b) An outer supersonic potential flow field 


Since an infinitesimal compression wave is assumed, the upstream 
pressure propagation is expected to perturb the boundary layer just 
slightly, thus imposing small disturbances upon the outer flow. This 
justifies the linearization of the supersonic outer flow. (Results of 
this theory will of course be applicable only for fairly weak shock 
waves.) The flow field has been assumed to be separated into two dis- 
tinct regions. In reality a continuous variation in velocity from zero 
at the boundary to some main-stream value a short distance from the wall 
is known to exist. The actual distance at which the velocity reaches its 
main-stream value is rather indefinite; hence the so-called boundary- 
layer displacement thickness 6*(x) is defined. Physically it represents 
how much the potential flow streamlines are deflected because of the 
reduction in mass flow caused by the retarded velocities near the wall. 
The viscous effects are thus characterized by a layer whose thickness 
is 6*(x), and the link connecting the outer flow to the boundary-layer 
flow is the condition that the line y = 8*(x) is a streamline of the 
outer flow. The treatment of the boundary layer in the present study is 
based upon a momentum-integral equation previously derived by Howarth 
(reference 13). The deflection of the boundary layer is expressed in 
terms of the vertical velocity component along the edge of the layer and 
this relation is used as a boundary condition for the outer flow. Hence 
the over-all problem is then reduced to the solution of a purely super- 
sonic potential-flow problen. 


LAMINAR -BOUNDARY-LAYER EQUATIONS 


The Prandtl boundary-layer equations for the two-dimensional steady 
flow of a viscous compressible fluid are: 
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Here u and v denote, respectively, the x and y components of the 
velocity and p, P, and T, the pressure, density, and temperature of 

the fluid. In order to simplify the problem, Howarth (reference 13) 
assumes that the coefficient of viscosity H is proportional to the 
temperature and the Prandtl number is unity. These assumptions, of 
course, introduce some approximation, but according to the results of 
Emmons and Brainerd (references 14 and 15), the boundary-layer character- 
istics depend very insensitively on the form of p(T) and the change of 
Prandtl number for moderate Mach numbers. Therefore, the simplifications . 
brought about by this assumption justify the slight loss in accuracy. 


By introduction of the stream function y 


dy 

pu = po >= 
° y 

= on 9Y 
pv Po SE 
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the momentum equation is transformed into 


EX 3x dx ax due E y 32x | - 33x 
(dy')3 


an (dyr)? d dx [Te eae? ^ (ay!) 


28e? 


where Po; Tos and v, stand for the density, temperature, and kinematic 


viscosity of a standard condition, the stagnation point, say, and are 
constants, and the subscript e indicates the condition at the edge of 


the boundary layer. In the case of a thermally insulated plate 


8 = 0, the condition of unity Prandtl number admits a special solu- 
W 


tion on the plate (reference 16) 


1 2 Ll 2 
Cpl +5 u = Cple + 5 Ue 


Upon elimination of T, there results the equation: 


32x dX 32x dX 33X 


— O nd 


axt Oy ày (y) Bet? (yi 


€ o NC ME e 
2a. — (3y!)* 


(1) 


d 
Thus, it is seen that when = = O, that is, at constant pressure, this 


equation is identical to the equation for the incompressible fluid. 


Let 5' be the boundary-layer thickness in the transformed plane, 
and define the corresponding displacement and momentum thicknesses by 


5!* Ti (7 = a dy! 
O RE 
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Then, by integrating eguation (1) across the boundary layer, the egua- 
tion relating the boundary-layer thicknesses with the velocity is 


+ gt + |l + 5 = (2) 
dx! "PL 2 2 u. dx! 2 1 
28e de e Le Oy W 


The similarity of this equation to that of the incompressible flow sug- 
gests the applicability of the Kårman-Pohlhausen procedure. 


If the velocity profile is given by (reference 13) 


u — 
ug = F(n) + AG() 
where 
F = 27 - 273 + në 
X 2 x 
G=zn(1-n) 
_ (51)2 y-1, AÀ du, 
a fi + 22 Me ax 
with 
1= E 
mo = 22 


it follows that 


b'* = (36 - A) (3) 


— — —- 
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Substitution of the quantities of equations (3) into equation (2) 
gives rise to a relation between the boundary-layer thickness 65' and 
the velocity uelx), with a parameter A. Now, in the present problem, 
as the flow in the boundery layer must be compatible with the potential 
field, the flow deflection d6*/dx of the boundary layer determines the 
flow in the outer field. Therefore, 5' should be eliminated. For this 
purpose, let the displacement thickness be introduced: 


(Cel E) (ao 


By means of the temperature-velocity relation and Pohlhausen's profile 
this will yield 


1/2 p A? 
5 (=~) T. É + 5 Me (a + t 3 2520 ^ * 252 36)" 


Now, from the second of equations (3), 5'* is already expressed in 
terms of 6' and A. The elimination of 6'*, then gives 


" 2j); à us 
Ee nes. 
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(ros me 112 r (2 ma. 
2 630 2 10 Po Te (4) 
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Theoretically, this equation defines 6' as a function of 6%. 
Only in the case of small pressure gradient, namely, 


Ue = VUl+u +...) 
ut << 1 


however, is a simple analytical solution possible. For by retaining 
only linear terms in u' and its derivatives there is obtained (refer- 
ence 17) 


2 
DL so 9| a ad O i563 data 
E 2520 © 5) = 315 6300 630 g^ u c 


+ 


2-7 1-37 
g2(7-1) _2-7 Mo20207-1) u!| 8% (5) 
2 : 
U y AE se 1 
where Mo =„ and o-l- 2 Mo . Moreover, since u' is assumed 


to be small, that is, for weak interaction, the boundary-layer thick- 
ness 6' must also differ from the unperturbed value by a small amount. 
Therefore, to the first-order approximation there results the solution: 


2-7 " 
o! = mee a 5% - 


315 6306 
o. 4-37 
263 (7 - 11M (63 _ 37y2 2-7 m.2o2(7-1) 263 — 30 asx + 
630 6-57 6300 315 2 6300 315 
(7-1) 
3(2- 
-2 ad -h t 

v; ns? . ES. - Any POD au (84) (6) 
Vo (3 2520 945/\6300 315 dx 


Let 5¿* be the unperturbed displacement thickness and write 


5% = 5,*(1 + A) 


a m e e P 
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Substituting equation (6) in the momentum equation (2), replacing x! 
by x, and keeping only first-order terms, one obtains, by taking into 


E 
n account 5¿% o = Constant, the following lihear equation (refer- 
x 
ence 17): um 
dA. uu du aut 0c 
od A a du! au! 
g, dê ka + 383 ge * 78h 3 (7) 
where 
= 2 263) Y lb Je 
* Re (37)(630) 263 
go = -81|1 a - DA = +35 
S é acto SLP 
263 
15934 74588 1 m 
"X Y)Mo > . 9731 360047 g 
i eo ] -= Th 0] 
263 
1 
— 3 =3 
= -R 07-12 (630 „sli 
gw e? ses 263 O 
with 
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The momentum integral thus relates the displacement-thickness perturba- 
tion A to the velocity perturbation and its derivatives along the edge 
of the boundary layer. 


SOLUTION OF BOUNDARY-VALUE PROBLEM OF OUTER FLOW 


Until now the study has been directed at obtaining an explicit rela- 
tionship for the boundary-layer displacement thickness as given in equa- 
tion (7). Now the extent of the interaction zone is only a fraction of 
the distance from the leading edge of the plate, hence 6,*(x) can be 


assumed to be constant over this range. Equation (7) then becomes a 
linear equation with constant coefficients. 


The problem is now one of determining the supersonic potential flow 
past a thin body whose thickness at any point corresponds to the perturbed 
boundary-layer thickness at that point. Consider the problem in the 
coordinates shown in figure 2. 


The undisturbed thickness is seen to correspond to y = 0. Strictly 
speaking, the undisturbed thickness should be at y = 55* in the physical 
plane. But since the outer flow does not depend upon the location of the 
origin of coordinates, the boundary conditions can be simplified slightly 
by choosing the axes as indicated in figure 2. Then, consistent with 
linear theory, the boundary conditions are satisfied at the undisturbed 
surface, namely, y = 0. 


The linearized differential equation for the perturbation velocity 
potential of the supersonic stream is 


2 2 
me oP _ 9D 0 (9) 
ot ane 


L and the general solution is given by 


where & 2 XL. and 7 = 
do 


Ø = ELE - mon) + g(& + mon) (10) 
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The incident waves will be taken as 


g(t + mm) = 0 when &£ + my <0 
(11) 
g(E + m,n) == 18 + 1) when E + mq >0 


This corresponds to a simple compression wave, of deflection angle €, 
incident upon the origin. (According to the above definition, the abso- 
lute value of e will be used in any computations.) 


Upstream Solution 


The incident wave causes the boundary layer to be perturbed upstream 
and downstream of the point of incidence. Since no disturbances can be 
propagated upstream into the supersonic flow, the physically possible 
solutions in region 1 will be waves of the form $4 = £(E - mon). The 
normal velocity component at the boundary must be zero, hence the condi- 
tion that the slope of the streamline is equal to the slope of the dis- 
placement boundary-layer thickness 5% is imposed. The velocity vector 


of the supersonic flow has' components [u(1 + u'),Uv'], where wu! = ag/at 
and v! - 36/31, so that the linearized boundary condition becomes 


aA _ op s 
a J at T O 
or 
da = «m f'(t) at 1 = 0 (12) 


Now since the boundary condition is on d A/dt rather than A, 
equation (7) is differentiated with respect to E, considering do* to 
be constant: 
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aA, dSĀ A deul, oa Aul . 
1 dA du! atu! 
Sunu ui 1 4 —— + emana | - 
1aAf au alu» ādu! 
81 ile ae * 983 o * 76h ua mn 


If the boundary condition (12) is substituted into equation (13), and it 
1 
is noted that at 7 =0, a f''(E), and so forth, an ordinary 


nonlinear differential equation for f(E) will be obtained. By line- 
arizing this equation: 


- yt (8) + gaf'''(6) + (82 + m,)£''(5) + gm ft (8) = 0 (1h) 


The coefficients in equation (1l) are constant for a specific case of 
Mach number and Reynolds number, hence the solution can be immediately 
written as 


TL E) = Ao T Ae 2* + Boe ls T coe 3 (15) T 
where the A's are roots of the equation 
gy? * 30° + (82 + Moo) A + gm, = O (16) 


For the cases under consideration equation (16) has three real roots, 
one positive and two negative. Let An >0, Aj < O, and Az <0. In 


order to have all disturbances vanish as E —+-o it is seen from 
eguation (15) that 


— nane porn ee = m ————— 000 A i — — e e mpfi ai — a M ci pi ma A — ———— 
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The upstream solution is then 


f(£.- mn) = m (17) 
or in the physical plane 
^o ( 
9 = pedof XY) 
| (18) 


X - Dy <0 


The solution in region 2 follows immediately. The eguation for the outer 
flow is linear, hence, by adding solutions, 


(19) 


Downstream Solution 


In region 3 the solution to the lineārized supersonic flow egua- 
tion is 


Pz = glt + mon) + h(E - mn) 


But g(t + Mon) is known since it is the incident wave, hence use of the 
boundary condition at the interface will result in an ordinary differen- 


tial equation for h(E). The nondimensionalized boundary condition takes 
the form 


ar = -mn'(&) - € &t n = O (20) 
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since 9%. -e - mh'(&) 
On , 
Substituting equation (20) into equation (13), 
gyhlY + gza'*! + (go + mo)h'' + gm! = ge (21) 


It can be seen that the complementary solution will be identical to 
eguation (15). Hence the complete solution is 


AoE OME Mat 


h(t) = Do + Be + Be + Ce - E (22) 


ER 
m, 
where again A2 >O, Al <0, and A3 <0. Choose Do = 0 without 

affecting any of the physical guantities such as velocity or pressure 
since h(t) merely represents a velocity potential. To eliminate the 


possibility of any velocity becoming infinite as ££ ——» œ put E=0. 
Then the solution becomes 


ME AE 
h(E) = Be + + Ce 7 - 6 
M(E-men) Nagl E-men) (23) 
h(E = mn) = Be * Lus E - -HE - mon) 
Then the perturbed velocity potential in region 3 is 
Al A 
— (x-m y) S (x-m y) 
TS E 
3 Mo Šo (24) 
X- My > 0 


Determination of Constants 


The problem now involves four constants A, B, C, and Dl with 
three parameters Mo c, and Re. Additional conditions must therefore 


be imposed in order to determine the constants as functions of the 


ithe constant D is a measure of the downstream boundary-layer 
thickness as shown in the following discussion. 


——— — -— — e a ee me er q... Tee O or 


NACA TN 2868 23 


parameters. It is noted that the outer flow field was divided into dis- 
tinct regions by the lines OS and OM (figure 2), so that the boundary 
layers in regions 1 and 3 were treated independently of each other. 

These boundary layers are therefore related at the point of incidence by 
the following conditions: 


(a) The boundary-layer displacement thickness must be continuous 
at x = 0 


(b) The wall pressure must be continuous at x = 0 


Xo 
(c) Discontinuities, if any, must satisfy lim K dx = 0 
xo O “Xo 
where K = O denotes the Karman momentum equation for compres- 
sible viscous fluids 


Condition (a).- Compute the boundary layers for regions 1 and 3 and 
then match the displacement thicknesses at the origin. One has as a 
boundary condition 


dA, ob 
—— = —— at = 0 i= lor 3 
af dn T f 
52% 
Where Aj = zu - l and O* is the total boundary-layer displacement 
O . 


thickness in regions l or 3. 


From equation (18) one obtains for region 1 


d Ay ^56 
aE = -mjoAe 


Not 
A, = -m Ae + Constant 


All disturbances vanish far upstream so that Aj(-») = O. But since 
No > O, the exponential term vanishes as 8 3 -a, hence the constant 
equals zero. Therefore 
ALE 
A = -mÃe 


or 
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For region 3 one obtains, according to equations (24), 


d A . 
=> = -mpre 5 mCr3e3 


E 
A3 = mæ - wee d + Constant 


„ Let A3( co) = D and since the exponentials vanish at positive infinity, 
the constant equals D. 


Therefore 
A 
Az = Be le - m Ce 35 + D 
or 
x x 
AL BF ^3 oF 
93% = 88% |1 - m,Be - me + D 


Thus it is seen that for the total downstream thickness 


ES a 53% = o*(1 + D) 


By putting 61* = 03% at x=0, 
A=B+C - = (25) 


Condition (b).- Since the validity of the usual boundary-layer 
assumptions in the interaction zone has been assumed in the case of weak 
shock interaction, the pressure at the wall, as a first approximation, 
equals the pressure at the interface. The pressure at the wall as well 
as at the interface is therefore required to be continuous. As a conse- 
quence, the pressure jump due to the incident wave must necessarily be 
neutralized by a reflected expansion wave. This deduction is actually 
confirmed by experiments. 


The pressure in each region consists of an undisturbed pressure Po 
plus a perturbed pressure p' due to the incident wave; P = pg * p'. 


Since the undisturbed pressure is the same throughout the flow field, 
only the perturbed pressures, which, by small perturbation theory, 


4G 
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are p! = -PxJu', will be considered. The velocity potentials dy, do, 


and $4 have previously been found, and, accordingly, the pressures are 
given by | 


À 
Bex xm) 
pj' = - JM PANDE (26) 
M | A 
EA xy) Beg x-my) gs 
pal = M Pol BNE + Cie E (27) 
The condition for contimuous pressure at the origin is 
P1' = p4' 
hence 
ec 
Año = BM + CM3 - = 28 
2 1 Si E (28) 


This condition shows that the discontinuity between regions 2 and 3 cor- 
responds, in linear theory, to an expansion of the flow from regions 2 
to 3 by a wave of the same magnitude as the incident wave. Now since A 


and X3 are negative, the downstream pressure along the interface y = 0 
is 


| 2 
YM 
P3' ae ud Pot2e) (29) 


as E —> œ. This asymptotic value of the pressure is exactly twice the 
incident pressure rise; therefore, as anticipated, the incident wave is 
reflected as a regular reflection from a solid boundary. 


Condition (c).- The Kármán momentum-integral equation for compres- 
sible fluids is 


— ——À —— A pra HEP ca Pala 
ja mae E = —————— 
na. ee me RA —— A. ARA mi lg: A eo mm — T — M —— ED pos RR e p e 
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Where Ty indicates the wall shear stress. With the ordinary definitions 
of the displacement and momentum thicknesses, the equation becomes 


du 
d e 
ās (peue?o) + Peuçd* — = Ty 


From the outer flow 


thus 


From the previous notation it is seen that 


— d 2 DT 
K = É (pouco - 8*p) tr Ty 


As the boundary-layer regions 1 and 3 are independently considered, 
the pressure gradient and velocity therein might be discontinuous at the 
origin. In order that the discontinuities be consistent with the dynami- 
cal equations, it is necessary that: 


Xo | 
d 2 do* 
lim ax Pete 0 - 5*p) + p = - Tyl dx = 0 | 


or 


26 Xg Xo ās* Xo 
lim (Pete - 2) P P ax dx - Ty dx | = O 
Xo —? 0 -Xo -Xo -Xo 


Now 5% is required to be continuous at the origin, p I discontinuous 
but finite, and the shear stress Ty, finite; then, in the limit, the 
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two integrals vanish and 


x 


O 
Lim (renee : p) = 0 
x9 -X) 
remains. Therefore, in the limit, 
P3u3 63 = py uy “0, = 5*(0) (P3 = P1) = 0 (30) 


The condition that p, = P3 at x = y = O has already been imposed 
and hence equation (30) reduces to 


3237601 = p3u3"83 at x = Y = 0 (31) 


This equation, linearized in u' and its derivatives, becomes 


duz' duy' i 
O - ST) + cafe! <) Ga 
Where 
-l 
1 1 -2 O 
2- = = 
o < Es 
ag =- 57 9 sd O 
.3T 263 y cd 37. 263) 
315  630c 315 6300 
2-7 
-1 
= gol XL. 283.) 2 21203 y «d. 3 
aa Gr 6300 TUO v 


The pressures were egual at the origin, so that, in linear theory, the 
streamwise perturbation velocities are also equal. Therefore uy" = uz' 


at x = y = 0. 


— A ——————— mm I q e < A Cre — ag 
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Now since a4 #0 and dp # O, for the range of M, and Re 
studied | 


dui : 


duz' 
mr m at x=y=0 (33) 


It is therefore seen that, in the linear theory, continuous pressure 
also implies continuous pressure gradient. Since 


du] ' A 

dx 5% at x=y=0 
duz' 1 
“a = gx ord at x=y=0 


there is obtained from eguation (33) 


Año? = BA? + CAZÍ (34) 


Conditions (a), (b), and (c) lead to equations (25), (28), and (34) 
for four constants A, B, C, and D. But previously it has been shown 
that 


lim 63% = 8,%(1 + D) 
E ce 


Thus it is seen that D merely determines the downstream boundary-layer 
thickness, whereas A, B, and C determine the local character of the 
perturbation. Therefore, if the downstream thickness can be estimated, 
there will be a determinate problem. Assume for the present at least, 
that D is determined. The constants A, B, and C can then be solved 
for from the following equations: s 


Aha = BM + Gha - E 
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The. solution is 


A -(A3 + X) EN MA E 
(^e - Sr - Y B (No - ^3) (^e - X) ^» 


= Da + 22) L Pe 


i (^ = Mla 7 Ay) =) i (a - LP - X) C E i 


. “(Ap + M) (- 29) : Me E 
(^3 = MIS = ^3 Mo (^3 E A1) (^e E A3) ( a.) 


Evaluation of D.- It is seen from equations (35) that the local 
perturbations are determined once the constant D is specified. As 
it has been shown that D characterizes the thickening of the boundary 
layer through & shock, it can be determined approximstely by the 
following consideration. 

When & weak shock is incident upon the boundary layer, it is &ssumed 
that the Kármán momentum-integral equation 


ie ld one y 2 + Xl Te NE EN 
dx Pa dx 6 Ju, dx | p 
ee 


is valid throughout the disturbed boundary-layer flow. For a given pres- 
gure distribution, the growth of the boundary layer is governed by this 
eguation. When a shock is incident upon the boundary layer, it has been 
theoretically predicted and experimentally verified that the flow up- 
stream of the point of incidence is separated over a considerable portion 
of the disturbed flow; hence the shear stress becomes relatively unimpor- 
tant. Downstream of the point of incidence, since the flow is dominated 
by pressure forces, whether the flow is laminar or turbulent the shear 

is known, from experimental results, to be very small (reference 18). 
Hence for the present approximation, the shear stress can be neglected. 


Furthermore, if the ratio 8*/@ =H is regarded as a parameter, 
the momentum equation gives the momentum thickness as a function of u(x). 


Now, it is known that H is substantially increased by the shock, and 
to a lesser extent by the compressibility. In going through a shock, H 
first increases because of the, shock and then decreases because of the 
drop in Mach mumber. Hence, as a first approximation, H can be taken 
to be some constant average value over the interaction range. (This is 


- = — 1 —— + rr arre MAA TA —— À AR ls B aro apne, Se < e Ge A eer ——— À ——— L mR 
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similar to the procedure mentioned, by Nitzberg and Crandall (refer- 
ence 19).) In particular, H will be taken as the value ahead of the 


interaction. 


The momentum-integral equation now becomes: 


or 


Integrating between regions 1 and 3, stations far upstream and downstream 
of the point of incidence, 


2+H 2+H 


OT 
2+H 
ae = PL kz (36) 


For laminar flow over a flat plate where p «T, Py = 1, 
and (OT/Oy), = O, the formula given by Lees (reference 20, p. 119) 


for H can be taken: 


H = 2.50 + 3.5022), 


The thickening predicted by equation (36) is shown in figure 3. 
Values of the density and velocity ratios were computed exactly and to 
first order in the deflection angle. For small deflection angles, the 
agreement between the "exact" thickness ratio and that obtained from 
linear theory is quite good. "The crosses on figure 3 are the experi- 
mental values for various Reynolds numbers obtained by Barry, Shapiro, 
and Neumann (reference 4), and there is fair agreement between theory 
and experiment. Actually, the experimental data are visual estimates of 
the ratio of boundary-layer thicknesses, so that for a comparison between 


— ~ ——— ms a — — —— n a, A e Ee GA A - ee ee e — 0. —— 
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theory and experiment, the ratio of boundary-layer thickness to displace- 
ment thickness would have to be assumed to be the same upstream as well 
as downstream of the point of incidence. In view of the assumption and 
the allowable experimental error in the visualization of schlieren photo- 
graphs, the experimental values must be regarded as qualitative. The 
theoretical values of the thickness ratio are near the lower limit of the 
experimental values. This is probably due to the fact that the value 

of H was underestimated, since an average H through the interaction 
range would be larger than the initial value. 


Evaluating D according to the linear theory, 


83% piu 
1* P343 


and — 
ey Mo, 
P3 Moo 
u 
Us Mo 
hence 
e+H 
2 = l+ (2 + E - Mo jos 
3143 Moo 
But 
52% 
x = 1+D 
1 
Therefore 
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32 
Now since 
H = 2.50 + 3.50(? 5 LIM? 
D = ES (4.50 = 0.3M) (37) 
Effects of Downstream Thickening 
Consider the constant A governing the upstream flow: 
oe -(a3 + X) (E) R A1A3 (-2) 
Pe - 23)(de - Ay) \ 9 Qe - 23)(de = XM) To 
or 
sec: "(13 + X) MA3 e : sd 
e) (ee MJA Be 


The coefficients are given in table I; hence it is seen that for all 
cases in the present range of Mach numbers and Reynolds numbers the con- 
tribution of D toward A amounts to less than lO percent. Thus it is 


seen that, in the present range, a precise knowledge of the downstream 
thickness is not essential to the determination of the upstream flow. 


The effect on the downstream flow is not quite so straightforward, 
since two constants B and C are involved. As an indication of the 
downstream flow, consider the pressure distribution along the wall. It 
has been seen that, far downstream, the pressure attains a constant 
value, namely the pressure that would have been anticipated had there 
been a regular reflection. Immediately behind the point of incidence the 
pressure is much lower than this end pressure since the flow has just 
undergone an expansion to the pressure that existed before the point of 
incidence. The manner in which the pressure proceeds from its value at 
the origin to its final downstream value is then of some importance. 

Does the pressure increase monotonically to its final value, or does the 
pressure at any point rise higher than the final value? Liepmann (refer- 
ence 3) has found experimentally that for the reflection of a shock wave 
from a laminar boundary layer there indeed exists a definite downstream 
overcompression. The pressure rises sharply to some value higher than 

l 


SS ——— i, rt: 9 ——— — —À —— m T ^ 


A ———— lições 


9Q 
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the end pressure and then tapers off. The experimental findings of Barry, 
Shapiro, and Neumann (reference 4) also confirm this observation, although 
their results show the overcompression to be slightly less pronounced. 


- Consider the behavior of the theoretical downstream pressure distri- 
bution. Along the interface the total pressure is given as 


x 
^ Bx ^ S x 
P3 = Po - MM 2po|Bme LO + Cage 379 - =: (38) 
and 
27M e 
BS PG 1 + Mo, ) (39) 
x 
as —— = E — > 00 
Šo 


so that far downstream, consistent with linear theory, the correct pres- 
sure behind the regular reflection is obtained. Then | . 


me hg + M) sh AM n 
K TEN 


wherein € > 0 and D>0. Now ha >0, XA1<0, and Az < 0, where Ao 
end A; are of the same order of magnitude, while ^3 is about ten 
times as small. Hence B <0. 


Consider 
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Now the second term is always negative in the Mach number range considered, 
whereas the first term may be either positive or negative depending upon 
the sign of ha + M. It turns out that for M, > 2.4 (approx.) the 


first term is positive, hence C > 0. For M, < 2.4 (approx.) the first 
term is negative, but 


Qo + X) : Mo .50 - 0.3M,? 


(s - a) (de - 73) (3 - aa) (Ne - Ag) is 


hence C is positive. Thus for the whole range of Mach numbers and 
Reynolds numbers, C > 0. 


The gualitative behavior of the downstream pressure distribution 
will now be established. The downstream pressure along the interface is 


X K 
2 AEK AB E 
Pz = Po + MP, = + YM, p,|-Bhje 190% - CAze 39o* 


The first two terms represent the constant end pressure while the third 
actually describes the variation of the pressure. Consider the curve 


x x 


A Bo A Pç 
Sex Sk 
y = -BMe-? - CMe 390 


It is known that B «0, C>0, A <0, and ^4 <0; hence -B <0 
and LA >0. The constants B and C are of the same order of magni- 
tude, but Aj is about ten times as large as 73; hence 


[an] > oal 


Therefore, when x = 0, y is negative and, when x becomes large, y is 
positive but small. The point at which the maximum occurs is given by 


2 

x) — ER E log BALD 

(55% a BN Cvm CA? 
max 
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If 


A mr pots 
y = -BAje 190% . Cge 390% 


is now interpreted as the pressure variation, it is immediately seen that 


the theoretical pressure distribution exhibits the overcompression that 
is observed experimentally. 


RESULTS AND DISCUSSION 


It has been seen that if the downstream thickening is estimated, 
the constants A, B, and C are completely determined as functions of 
the flow parameters. With the outer flow thus determined, the boundary- 
layer growth and pressure distributions can now be computed, and hence 
the effects of Mach number, Reynolds number, and shock strength upon the 
upstream influence and the location of the separation point can be 
studied. (Of course, when the results are compared with experiment, the 
Reynolds numbers must be low enough so that the boundary layer will 
remain laminar, prior to the interaction, in the corresponding experi- 
mental case.) 


Upstream Influence 
In order to estimate the upstream influence of the interaction, 
define a length in which the pressure on the boundary-layer displacement 


thickness decays to a specified fraction of its amplitude at the origin. 
This length would then be determined solely by the exponent Ao. For 


example, consider the upstream boundary-layer disturbance: 


A 
81% = -m_AS te 200% 
1* = -mA5g*e 


At the origin 


5,*(0) = -m,A5,* 


= —- mm r e A a ee c 
— ee N 
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If xa/5o* is defined as the distance required in order that 
6]* = a( -m,A89*) at xa/9o* where a < 1, 


x 
Lear's 
a (-moABo*) = -mot te 200 
Therefore 


Os” ^o e CL 


Since a is a constant it is seen that the upstream influence is then 
inversely proportional to Ao. Values of ha have been plotted in fig- 


ure 4. (Values of ^, and 24 are plotted in figures 5 and 6, respec- 


tively.) For fixed Mach number, xa/30* increases with increasing 


Reynolds number. For large Reynolds number, this dependence can be 
deduced from the equation: 


S Bats (Be Bh ea 
gh gh Eh 8) 


The coefficients are functions of Mach number and for the range of Mach 
numbers considered, the roots are real. When the Reynolds number is 
large, it can be shown that 


ho MN 
-8Ļ 


where gy (Re, M.) <0 and gy «Re. This is in agreement with Lees' 


result. Consequently, for large Reynolds numbers, xa/8,* « Rel/2, For 
fixed Reynolds numbers, the upstream influence decreases with increasing 
Mach number as indicated in figure 4(b). If the disturbance is considered 
to decay to, say, > percent of its amplitude at the origin, it is found 


that, at the high Reynolds mumbers, and Mach numbers about 2 or less, 
the upstream influence is of the order of 30 boundary-layer displacement 


thicknesses. 
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The result that the upstream influence increases with Reynolds num- 
ber may &ppear to be disconcerting since the viscous effects are expected 
to be more prevalent &t the lower heynolds numbers and hence to produce 
a greater influence upon the upstream flow. This &pparent paradox arises 
because the upstream influence has been measured in multiples of & length 
which is also dependent upon Reynolds number. If the absolute values of 
upstream influence are considered, 


fà pel/2 
do? ( 
But 
hence 
Xa * 5,* Rel/* 
Or 
1/2 


Xa a Re. 


Thus, measured on an absolute scale, the upstream influence increases 
with decreasing Reynolds number. 


As defined above, the distance of the upstream influence is depend- 
ent solely upon Mach number and Reynolds number. It can be redefined so 
that it will also depend upon the shock strength. This is accomplished 
by defining the upstream influence to be that distance at which the dis- 
turbance decays to a given fraction of the undistrubed value. The per- 
turbed boundary-layer displacement thickness ahead of the point of inci- 
dence is given by: 


2 x 
= 209% 
07% = -m_A5,*e 


where 6 ,* is the undisturbed displacement thickness. The upstream 
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influence is then the distance xa/do* to the point where 87% = bå. % 
and b < 1: 


xX 
x 
b85* = -m,AB Xe 250 
Kā 
Ao Ex = loge [s "wd 
O 
šā iu mn 10 we 
KAARE “8 b 


Now A is negative and decreases with increase of the deflection angle, 
so that for fixed Mach number and Reynolds number, the upstream influ- 
ence decreases with shock strength. If the upstream influence is mea- 
sured to the point at which the disturbance decays to 5 percent of its 
undisturbed value, b = 0.05. In figure T, xa/do is plotted against 
deflection angle. For ¢ > 29, the results seem to agree fairly well 
with the experimental values presented in figure 10, reference 2. For 
larger deflection angles, however, the theoretical values are too low in 
comparison with the experimental values. This might be due to the fact 
that the linear theory becomes less accurate as the deflection angle 
increases. 


* Boundary-Layer Separation 


Since the outer potential flow is known, the point of separation 
can be computed on the basis that separation occurs when (du/dy), = O 


Now ou/dy is proportional to Ou/y' so that one can just as well use 
(du/dy'),, = O as a criterion for separation. In view of the assumption 


of a Pohlhausen velocity profile, in the transformed plane 


du Ue A 
I s = lo EX 
A sr ( t ē) 


hence, as in the incompressible case, separation occurs at A = -12. 

Now it is known from experience in the incompressible case that, in 
regions of retarded flow, the Kármán-Pohlhausen method gives values of 
the skin friction that are too high and consequently predicts separation 
too late or not at all. This feature of the method is also to be 


— -— vii ar Fa eS 


NACA TN 2868 39 


expected in the compressible case, since it was shown by Howarth that 

the effect of compressibility is eguivalent to exaggerating the pressure 
gradient in the incompressible case (reference 13). In particular, 
Stewartson has used the Kārmān-Pohlhausen method to treat the case of 
flow against a linear pressure gradient and has shown that the predicted 
distance from the leading edge to the separation point is an overestimate 
of the actual value (reference 21). 


In the present problem, where the flow upstream of the point of 
incidence is subject to a positive pressure gradient, the same overesti- 
mation of the separation point will, of course, be expected. But the 
effects of Mach number, Reynolds number, and shock strength upon the 
separation point are of primary interest, so that inaccuracy in the abso- 
lute values is unimportant. 


It is known from the boundary-layer theory that the parameter that 
defines the separation point is 


~ 


2 
NC OMNE ATTE 
A = Va dx l+ 2 Me 


By neglecting products of u' and its derivatives, 


A= [1573/46 (0/2) aut /ax) 


Now, from equation (6), to the order of approximation, 


3-27 
Re Šo* y-1 du! 
==——————— Y mc 
E = 
315 6300 


Along the edge of the boundary layer, from the upstream solution (see 
equation (18)) 


2 x 
aut a 2 MBF 
dx “5 


E € ==. -= — —€———————Á—— P A m mm a 
o o Loc 7 E ui 
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so that 


E 
Ne ac MBF o* 


ESAE 
315 6300 
For separation 
3-27. *s 
2 A . m 
-12 si E L 
31 263 Y 
Í 315 6300 
Xg 1 A Re hp? 
id loge (Tip B (40) 
where 
3-27 
pl 
B = 


= 2 
37 963 174 
315 6300 
and Ao is positive, hence separation ahead of the point of incidence 


will be predicted only when - = A Re Mo B > 1. 


Contrary to the separation phenomenon in ordinary boundary-layer 
flow, the location of the separation point is influenced by the Reynolds 
number. This is easily understood when one considers the fact that 
separation of the flow is controlled by the pressure gradient of the 
outer flow; according to boundary-layer theory, the outer flow depends 
solely upon the geometry of the body. Consequently, the location of the 
separation point is independent of Reynolds number. In the present 
problem, however, an outer flow that is compatible with the boundary- 
layer flow has to be found. This relationship is expressed by the condi- 
tion which requires that the direction of the potential flow be the same 
as the slope of the displacement thickness. Since the boundary-layer 
growth depends upon the Reynolds number, the outer flow, and eventually 
the separation point, mist also vary with Reynolds number. 


6Q 
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In order to determine the effect of shock strength upon the position 
of the separation point, the constant A must be examined: 


A =-28 a MM : HS 4.50 - > 
Moo (^e = A3) (X2 = Xi) Mo ~ A3) No 7 M1) Tho 


For the range of Mach numbers and Reynolds numbers considered, the terms 
in the bracket are positive, since ^54 > O, A < O, and Az < O. There- 


fore as € is increased A becomes more negative; consequently 


- 5 A Re ^o*8 becomes more positive and hence the separation point moves 


upstream. Thus an increase of shock strength increases the distance 
between the separation point and the point of incidence. 


To express xs/d0* explicitly in terms of Mach number and Reynolds 


number would be difficult; hence this relationship will be presented 
numerically by varying, separately, the Mach number and Reynolds number. 
In figure 8, lines of constant Reynolds number are plotted in the 
xg/50"5€ plane, and the separation point, measured in multiples 


of 59%, moves upstream with increasing Reynolds number. This is not 


too surprising since the same behavior occurs with the upstream influ- 
ence. As the upstream influence increases, the "self-induced" pressure 
gradient will begin farther upstream and hence separation will occur 
farther upstream. For fixed Reynolds number, the upstream influence 
increased with decreasing Mach number, and accordingly (see fig. 9) 

the separation point, measured in multiples of 5¿%*, moves upstream with 


decreasing Mach number. 


Before closing the discussion on separation, the importance of D 
in the determination of the separation point should be discussed. The 
coefficients in the expression for A are such that, for fixed Reynolds 
number, Mach number, and deflection angle, the magnitude of A increases 
as D increases. The separation point has been seen to move upstream 
as A is increased in magnitude, hence an overestimation of the down- 
stream thickening would result in a slight overestimation of the upstream 
distance to the separation point. In figure 9 the variation of separa- 
tion point for different Mach numbers has been plotted for fixed Reynolds 
number. The location of the separation point for the case D = O has 
also been plotted in this figure. These curves thus give the greatest 
lower bound of the separation distance since it is known that actually 
D >0. A comparison of results reveals that D has a very small effect 
upon the location of the separation point. For the higher Mach numbers 
in the range investigated, the percentage difference between the two cases 
may be fairly large, but the fact that the separation distance is mea- 
sured in multiples of a boundary-layer displacement thickness must be 
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considered. In an actual experimental measurement of the absolute dis- 
tance to the separation point, this difference will be extremely small. 
Hence, for practical purposes, one can consider D = O when estimating 
the separation point. The results obtained under such conditions will 
then yield a slight underestimation of the separation point. 


Pressure Distribution 


The pressure disturbance along the wall (fig. 10) decays exponen- 
tially from a definite value at the point of incidence to zero far 
upstream of that point. Downstream of the point of incidence the pres- 
sure rises to a maximum value before dropping to the value corresponding 
to regular reflection. This downstream overcompression has been observed 
experimentally, and it appears to be a characteristic feature of shock- 
wave interaction with a laminar boundary layer. The Lees' theory, as 
mentioned previously, failed to predict this downstream behavior. This 
is due to the fact that an incomplete solution for the pressure was used 
in the determination of the pressure distribution. The boundary layer 
was divided, longitudinally, into four regions, and it was assumed that 
the solutions to a third-order differential equation were valid in each 
region. In the two regions that extended to positive and negative infin- 
ity, certain solutions could be rejected since they became infinite at 
the ends of their respective regions. In the finite regions, however, 
the complete solutions must be retained. The incompleteness of the Lees' 
theory, then, arises from the fact that only one term of the general 
solution was used in each of the finite regions. 


The linear theory yields pressure distributions that are similar 
irrespective ‘of the size of the deflection angle. Now for very small 
angles, the experimental results exhibit the general behavior predicted 
by theory. In figure 10 the experimental values of the pressure distri- 
bution have been plotted for e= 1° and M = 2.05. (In this case, 


separation has probably not occurred since the wave is quite weak.) The 
values of the pressure ratio were taken from figure 13 of reference Ļ 
and converted to the scale indicated. The upstream portion of this 
curve can be well represented by an exponential curve, thus verifying, 
at least for this case, the predicted exponential pressure rise. For 
larger deflection angles, the experimental pressure distributions are , 
characterized by the familiar pressure "bump" ahead of the point of inci- 
dence, thus indicating that separation has occurred. The linear theory 
is inadequate, as regards predicting this upstream behavior; hence for 
large deflection angles the theory mst be modified, possibly by taking 
account of the second-order effects upstream of the point of incidence. 
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SHOCK-WAVE INTERACTION WITH A TURBULENT BOUNDARY LAYER 


Method of Approximations 


When the flow in the boundary layer is turbulent, the Kármán 
momentum-integral eguation will take on the same form as for the laminar 
boundary layer except that the guantities involved are not the exact, but 
average, values. However, since general relationships betveen the shear 
stress and the mean velocity in turbulent layers have not yet been estab- 
lished, rigorous treatment of the problem at this time is, of course, not 
possible. Here, for the purpose of exhibiting the characteristic dif- 
ference between the laminar and turbulent cases, certain approximations 
are proposed. 


The momentum-integral equation expresses the shear stress at the 
wall in terms of the growth of the boundary-layer momentum thickness and 
the velocity gradient. Using experimental results as a guide, the rela- 
tive importance of these terms can be estimated. Experimental results 
of Fage and Sargent (reference 18) show that the shear stress in front 
of the shock is practically constant; behind the shock, it is very small. 
Therefore, unlike the laminar case, the shear stress to the first appproxi- 
mation, both in front of and behind the shock, can be regarded as con- 
stant and hence has no effect on the perturbed flow. It follows then 
that, in the case of a turbulent boundary layer, the growth of the momen- 
tum thickness is influenced pM M pressure changes due to the 
presence of the, shock. 


To simplify the problem further, it is noted, for practical pur- 
poses, that the shape parameter H =5*/0 is relatively insensitive to 
change even though there may be a considerable adverse pressure gradient. 
In the transonic case, where there is a normal shock in the local super- 
sonic region and a large change in H is anticipated, H at its maximum 
is only increased by a factor of about 1.2 (see reference 2). The reason 
for this is possibly the fact that the increase by the shock may be 
counterbalanced by a decrease due to compressibility effect. Moreover, 
in the momentum-integral equation, the coefficient of the velocity gra- 
dient is positive and usually greater than unity within the Mach number 
range; hence one can ignore the slight variation in H and consider it 
a constant. By the same token, the Mach number in this coefficient can 
also be considered to be constant. (The Mach number enters when the 
density is eliminated in terms of the velocity.) Consequently, the 
momentum-integral equation reduces to: 


t 
dA, a3 du -0 
dx dx 
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where 


az =2 +H- Moe 


Solution for Outer Flow 


Similar to the previous case, for the upstream region _ 
$1 = f(E - mm). By applying the original condition dA/āt = dø/dn 
at my = O, there is obtained from the linear theory 


f= cQ aL E-men) + Constant 


Since m,/a3 > 0, the disturbances vanish far upstream. When a3 is 


evaluated the exponent Med/ 3 is found to be about one order larger 
than As of the laminar case. This is the well-known experimental 


result that there is very little upstream influence in the turbulent 
case. This, therefore, confirms the hypotheses made in the previous 
section. 


For the downstream solution 


$3 = (E - mn) - (E + men) 


The original differential equation remains the same; hence the down- 
stream solution would be 


da = d og Emon) | = t + Constant 


It shows that if the solution 1 is contimed to the downstream side 


the velocities would become infinite at positive infinity and must be 
rejected. There is therefore a principal difference between the laminar 
and turbulent cases. In the laminar case several solutions for the outer 
Plow were obtained and the solution appropriate for either upstream or 
downstream could be chosen. In the turbulent case there is no choice 
since there is only one solution for the outer flow, which, if continued, 
fails at positive infinity. This indicates that linearization of the 
flow is incapable of accounting for the flow near the point of incidence 
where the nonlinear effects become important. Since large changes in 
flow velocities can be brought about only through a shock, in the pres- 
ent problem a reflected shock must be considered. 
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From general considerations it can be argued that if one puts a 
shock wave between regions 2 and 3 and considers perturbations of a regu- 
lar reflection it may be possible to obtain a solution which will satisfy 
the conditions at infinity. By perturbing the regular reflection, there 
would be, in each region, undisturbed quantities plus their perturba- 
tions. Consider the total potential in region 3: 


$3 =E + 2(6 - my) + g(& + my) 


At infinity f' = g' = 0; along the reflected wave E - man = O, to first 
order, the velocities are constant; hence 


f' + g' = Constant 


But to first order, along the wave, f'(0) = Constant; hence 
g'(£) = Constant. Therefore, g'(&) is constant throughout region 3. 
But g' = 0 at infinity; hence g' = 0 in region 3. 


Then at the boundary layer à A/āt = -m3f'(§) where f is nondi- 


mensional, and the substitution into the momentum-integral equation 
yields 


-m4f' + azf!! = O 
ma la 
f = Constant + foe 3/ 38 
m GE 
i= fo = al 3 
ES 
Now f' =0 at infinity; hence one must choose fo = O, and thus 


f = Constant 


This means that behind the reflected shock, the flow, to the first order, 
is uniform. Since the flow behind the incident shock has undergone com- 
pression up to the refiected wave, it can be seen that, to satisfy the 
condition at infinity, the reflected shock must be followed by a very 
rapid expansion. Otherwise, the pressure after the second shock would 
have been higher than that after the regular reflection. It is there- 
fore concluded that in the case of a turbulent boundary layer the inci- 
dent shock is reflected as a shock and behind the reflected shock there 
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must be an abrupt expansion so as to cancel the overcompression brought 
about by the train of upstream compression waves. The existence of the 
sharp expansion behind the reflected shock was also confirmed by experi- 
ments (reference 9). The pressure distribution and boundary-layer 
growth, to the first order, appear to be discontinuous as shown in fig- 
ure 12. 


Downstream Thickening for Turbulent Boundary Layers 


The downstream thickening can be estimated by use of the formula 


e+H 
53% P1 uz 
97 * E Paus 


Now the difference between the laminar and turbulent cases is that H 
differs for the two cases. Since, in general, Hy, < Hz, the downstream 


thickening in the turbulent case is less than the thickening in the 
corresponding laminar case. Empirical relations mist be relied on to 
estimate H because of the lack of knowledge of the compressible turbu- 
lent boundary layer. According to Nitzberg and Crandall (reference 19), 
for local Mach numbers greater than the free-stream Mach number, but 
less than 1.4, the compressibility effect is well approximated by 


s-a EIL) 


Use this relationship to estimate the downstream thickening and 
then compare the results with the existing experimental data. By 
assuming a 1/7-power law for the incompressible profile, 


H = 1.29(1 + 0.1M?) 


The predicted thickening is shown in figure 13. It shows that the down- 
stream thickening increases as the Mach number is decreased. This 
behavior is also present in the laminar case. 


Using this value of H, it was found that the predicted thickening 
for My = 1.44 and ¢=4.5° is larger than the visually estimated 


thickening in the experimental case. On the other hand, the predicted 
thickening for M,=2 and ¢=6 is slightly less than the visually 


estimated thickening in the experimental case. It should be noted, how- 
ever, that, for M > 1.k, a relationship has been uséd that is supposedly 
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valid for M< 1.4. Moreover, this relationship merely accounts for the 
effects of compressibility. It would appear that some effect of the 
shock must also be included in the determination of H. However, at 
this time, owing to the lack of exact information, it must be neglected. 


Consequently, it appears that if the relation H = 1.29 (1 + o.) is 


used for Mach numbers not too much larger .than 1.4, one can expect to 
have & minimum estimate of the downstream thickening. 


SUMMARY OF RESULTS 


An investigation of the reflection of a weak shock wave from a 
boundary layer along a flat plate yielded the following results. 


Laminar Case 
l. In all cases investigated, the pressure along the wall overcom- 
presses downstream of the point of incidence. The pressure disturbance 
decays exponentially from a definite value at the point of incidence to 
zero far upstream of the point of incidence. Downstream, the pressure 
rises to a maximum value and then falls, gradually, to the constant 
value corresponding to regular reflection. 


The exponential pressure rise appears to be verified in the case of 
a shock deflection angle of 3% since separation has probably not 
occurred. For larger deflection angles (the next larger angle for which 
there are experimental pressure distributions is 3°), the experimental 
pressure distributions exhibit the familiar pressure bump between the 
separation point and the origin. For these angles, a true comparison 
between experimental results and theoretical results cannot be made 
since the present theory does not account for the effects of separation. 


2. If the upstream influence is considered to be the distance to 
the point at which the disturbance has decayed to a specified fraction 
of its amplitude at the origin, the upstream influence, when measured in 
multiples of the boundary-layer displacement thickness 50%, is found to 


increase with increasing Reynolds number:  x/89* « (Re/m,)+/2 where x 
is the coordinate parallel to the flow direction, Re is the Reynolds 


number, M, = Ee l, and Mo is the free-stream Mach number. For 


decreasing Mach number, the upstream influence also increases. If the 
disturbance is considered to decay to, say, 5 percent of its amplitude 
at the origin, the upstream influence for Mo %2 and Re x 1500 is of 


the order of 30 boundary-layer displacement thicknesses. 


er aeo ra e-r a ——— mm 
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3. The pressure gradient is such that the boundary layer may separate 
ahead of the point of incidence. In the event of separation, an increase 
of shock strength, for fixed Mach numbers and Reynolds number, increases 
the distance between the separation point and the point of incidence. 

For fixed Mach number and shock-deflection angle, the separation point 
measured in multiples of Š9* moves upstream with increasing Reynolds 


number. For the case of M, =,1.44 and Re = 2000, the boundary layer 
separates for all flow-deflection angles ¢ > 1.03%. 


For fixed Reynolds number and shock-deflection angle, the separa- 
tion point, measured in multiples of do*, moves upstream with decreasing 


Mach number. 


For a complete determination of the constants of integration, an 
estimate of the downstream boundary-layer thickness was required. For 
the cases investigated, the effect of the downstream thickening on the 
outer flow is rather small. In fact, for practical purposes, this effect 
of downstream thickening can be neglected in the determination of the 
Separation point. The distance between the origin and the separation 
point will then be slightly underestimated. 


4, The present theory is applicable only in the case when a weak 
shock is incident upon a laminar boundary layer. In addition to the 
outer flow field, the boundary-layer displacement thickness has also 
been linearized. The latter linearization enables a linear differential 
equation with constant coefficients to be obtained for the perturbation 
velocity potential. This equation immediately yields the general solu- 
tion of the outer flow. Upon investigating the size of the perturba- 
tions, the maximum velocity perturbations are found to be about 10 per- 
cent or less. The thickness perturbations, on the other hand, are much 
larger, being about 30 percent near the origin; thus the linearization 
of the displacement thickness becomes questionable as the point of inci- 
dence is approached. Had the displacement thickness not been linearized, 
a rather complicated nonlinear differential equation would have been 
obtained from the boundary condition for the perturbation velocity poten- 
tial. A solution of this equation would be expected to yield more accu- 
rate results. However, it is problematical as to whether the equation 
could be solved without the imposition of additional assumptions which, 
in themselves, might nullify.any accuracy that the nonlinear boundary 
condition may provide. 


Turbulent Case 
In the case of shock-wave interaction with a turbulent boundary 


layer the upstream influence is found to be considerably less than in 
the laminar case. In addition, to first order, the incident wave must 


TQ 
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be reflected as a compression wave followed immediately by an expansion 
wave, so that the end pressure condition is satisfied. 


Cornell University 
Ithaca, N. Y., January 11, 1952 
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TABLE I 


COEFFICIENTS USED IN DETERMINING CONSTANTE OF INTEGRATION 
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Figure 1.- Simplified model of shock-wave boundary-layer interaction. 
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Figure 2.- Shock-wave boundary-layer interaction showing choice of axes. 
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(a) Various values of M, 


Figure 4.- Measure of upstream influence. xg/do* = - > loge (1/0); 
disturbance at xg/do* is a percent of disturbance at origin. 
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(b) Various values of R,. 
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Figure 4.- Concluded. 
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Figure 5.- Reynolds number against Al. 
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Figure 6.- Reynolds number against | 13. 
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Figure 7.- Upstream influence against deflection angle. M, = 2; 
Re = U8g*/Vg = 2000; and Rex = UxÃo * 51,600. 
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Figure 8.- Separation point against deflection angle. M, = 1.44; 
Re = U8 * / Vo. 
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Figure 9.- Separation point against deflection angle. 
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Figure 10.- Wall pressure distribution. M, 2; Re = USp*/V,. 
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Figure 11.- Perturbed boundary-layer displacement thickness. 
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Figure 12.- Downstream thickening turbulent boundary layer. 
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Figure 13.- Wall pressure distribution. 
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